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Abstract 

The production of the octet of baryons and mesons in e + e~ collisions is analysed, 
based on considerations of SU(3) symmetry and a simple model for SU(3) symmetry 
breaking in fragmentation functions. All fragmentation functions, D^(x, Q 2 ), describ- 
ing the fragmentation of quarks into a member of the baryon octet (and similarly for 
fragmentation into members of the meson octet) are expressed in terms of three SU(3) 
symmetric functions, a(x,Q 2 ), /3(x,Q 2 ), and j(x, Q 2 ). With the introduction of an 
SU(3) breaking parameter, A, the model is successful in describing hadroproduction 
data at the Z pole. The fragmentation functions are then evolved using leading order 
evolution equations and good fits to currently available data at 34 GeV and at 161 
GeV are obtained. 

PACS Nos: 13.65. +i, 13.87.Fh, 13.85.Ni 



1 Introduction 



The formation of hadrons from fragmentation of partons is of considerable current 
interest Jl|, [|. While the production of partons in any process — be it Deep Inelastic 
Scattering or e + e~ annihilation — can be calculated, perturbative QCD can only pre- 
dict the scale (Q 2 ) dependence of the process of fragmentation of quarks and gluons 
into hadrons; the fragmentation functions themselves are not perturbatively calcula- 
ble and can only be modelled. Various models exist which attempt to explain the 
process of fragmentation ^ |], |5|]. Many computer simulations || [7| also exist and 
are in popular use. However, the role of strangeness suppression as well as isospin in 
hadroproduction is not yet clearly established. A clean channel to study such phe- 
nomena is provided by e + e~ annihilation experiments due to the fact that the initial 
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interacting vertex is purely electromagnetic in nature. These experiments have been 
performed at different energies || || |i~0f . 

We propose a simple model for a light quark (u, d, s) to fragment into an octet 
baryon or a pseudoscalar meson, using SU(3) symmetry of quarks and octet hadrons. 
All fragmentation functions are described in terms of three SU(3) symmetric functions 
a(x, Q 2 ), (3(x, Q 2 ) and 7(2;, Q 2 ) (one set for baryons and another set for mesons) and 
an SU(3) breaking parameter A which have been determined by comparison with 
data. The model is able to predict the ^-dependence of all octet baryons and mesons. 
There is good agreement with data at two different energies (corresponding to Z° 
and photon exchange) over most of the x range of available data. Hence, the overall 
success of the model does seem to indicate the existence of an underlying SU(3) 
symmetry between members of a hadron octet. 

The paper is organised as follows: in the next section, we develop the model for 
quark fragmentation into octet baryons and mesons. In Section 3, we fix our model 
parameters using data on some hadrons at the Z° pole and use the resulting fits to 
predict production rates for other hadrons. We find good agreement with data. In 
Section 4, therefore, we use leading order evolution equations to simultaneously fit 
data at two different energies corresponding to hadroproduction via Z° and photon 



exchange. We use both the standard DGLAP evolution equations |TT| as well as 
the modified leading log approximation (MLLA) [12|] to evolve the fragmentation 
functions to different energies and compare the model with data. Section 5 contains 
discussions on our results and concludes the paper. Some details for the expressions 
of the different quark fragmentation functions in terms of a, (3 and 7 are given in 
Appendix A. 



2 Cross Section and Kinematics 

We consider the production of hadrons in e + e~ annihilations via 7 and Z exchange. 
To leading order, the cross-section for producing a hadron h can be expressed [[T<| in 
terms of the unknown fragmentation functions, D^(xe,Q), as 

^c q D h q (x E ,Q) 
1 da q 



crtot dx E 
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Here c q are the charge factors associated with a quark qi of flavour i and can be 
expressed jnj in terms of the electromagnetic charge, e,, and the vector and axial 
vector electroweak couplings, v j = T 3i — 2e, sin 2 9 W and cij = T 3i as 



c = c v + a 

V 47TO 2 2 ( 2 , 
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In Eq. |I|, a sum over quarks as well as anti-quarks is implied. Here xe is the energy 
fraction, xe = -Ehadron/-Ebeam = 2E h / \fs. We shall also use the momentum variable, 
x p = -Phadron/ -Pbcam = 2Pft/ y/s] x 2 E — x 2 + 4m\/ s, where nih is the mass of the hadron 
h and Q is the energy scale of the interaction and is equal to \/s. We shall normally 
use x to mean xe unless otherwise specified. 

The fragmentation function, D^(xe, Q), associated with the quark q is the proba- 
bility for a quark q to hadronise to a hadron h carrying a fraction xe of the energy of 
the fragmenting quark. This is not perturbatively calculable from theory, although 
the scale (Q) dependence of these functions is given by QCD. Data from different 
experiments at different energies from yfs = 10-91.2 GeV exists on p, A, E and S 
octet baryon production as well as on tc, K, rj octet meson production 0, |||, [Mj. All 
available data measures the production rate of hadron plus antihadron. Due to the 
symmetric nature of the process e + e~ — ► qq, the resulting hadron and antihadron 
yields are equal. We therefore present results for the sum of hadron and antihadron 
yields in what follows. 

We can also re-express the cross section in terms of the octet and singlet fragmen- 
tation function combinations, as 

1 da h _ a E h (xE, t) + a 3 D 3 (x E , t) + a s Ds(x E , t) 

cr t ot dx E C Q 

i 

where S, D% and D$ refer to the singlet, and the two octet ((u — d) and (u + d — 2s)) 
combinations respectively, with a = (c u + Q + c s )/3; a 3 = (c u — c d )/2 and a 8 = 
{c u + c d - 2c s )/6. 

We now present our model for quark fragmentation functions. 



Pi( s ) = 

P2(S) = 



3 The Model 

We study semi-inclusive production of the octet baryons and the pseudo-scalar octet 
mesons in e + e~ annihilation process using SU(3) symmetry of the quarks and of the 
hadrons in their respective octets. The production of the entire meson and baryon 
octet is described in terms of SU(3) symmetric quantities. We study only light quark 
(u, d, s) fragmentation where the fragmenting quark (g^) is a member of the quark 
triplet (qi = u, q 2 = d, q 3 = s) and the hadron under study (hf) is a member of the 
baryon (or meson) octet (see Table |l|): 

q -> h + X , 

with X being a triplet, antisixplet or fifteenplet. Since the final states in these three 
processes are distinct, we can express all the quark fragmentation functions for the 



3 



hadron h in terms of the independent fragmentation probabilities a(x,Q), ft(x,Q), 
and 7(x, Q) where a (ft, 7) is the fragmentation probability when X = 3(6, 15). The 
corresponding probabilities for antiquark fragmentation are a, ft and 7 (in this case, 
there is an antitriplet fragmenting into an octet and X, with the X being a 3,6 or 15). 
These probabilities are also functions of xe- The quark and antiquark fragmentation 
into h\, i,j = 1, . . . , 8, in terms of the SU(3) probability functions, a, ft, and 7 are 
given in Table ^| Note that the functions a, ft and 7 for the baryon and meson octets 
are unrelated; they just correspond to the same underlying symmetry. 

Since the (more) massive strange quark is known to break SU(3) symmetry, we 
introduce symmetry breaking effects as follows: the fragmentation function is sup- 
pressed by an x^-independent factor A whenever a strange quark belonging to the 
valence of the hadron is produced. This means that all non-strange fragmentation 
functions of strange hadrons are suppressed by A. For example, or are sup- 
pressed by a factor A compared to Df or D^. Note that all (strange and nonstrange) 
sea fragmentation functions corresponding to a given hadron come with the same 
factor of A. 

There are 3 quark and 3 antiquark fragmentation functions for each hadron. Since 
there are eight baryons (or mesons) in the octet, this corresponds to a total of forty 
eight unknown functions for a given octet, which, in principle, need to be fitted to 
data. However, in our model, all the hadron fragmentations for a given octet can 
be described in terms of 6 functions alone, along with an SU(3) breaking parameter, 
A, thus leading to an enormous simplification in the analysis as well as dramatically 
increasing the predictive power of the model. 

We now go on to detail the model, first for the case of octet baryons and then for 
the octet mesons. 

4 Comparison with data at the Z pole 

We relate the model parameters and try to find a, ft, 7 (and the corresponding 
antiquark functions) as well as the SU(3) breaking parameter A by comparing with 
data on e + e~ — > Z° — > qq obtained at LEP || || . 

4.1 Baryon Fragmentation 

First of all, we observe that our model predicts equal rates of production of S° and 
(S + + £~)/2 (See Table ||) due to isospin symmetry. This is borne out by data ||; for 
instance, the multiplicities of (£+ + £~)/2 and S° at Q = 90 GeV are 0.091 ± 0.019 
and 0.071 ± 0.018 respectively and are compatible 1 within la. 

To obtain predictions for the other baryons, we separate the quark fragmentation 
functions into valence and sea parts by defining, as usual, 

a = a v + a s ; ft = fty + fts \ 7 = Iv + Is ! 

1 Our model does not account for isospin breaking effects which are small compared to SU(3) 
breaking effects but are known to exist; hence we are looking for agreement only to within this 
approximation. 
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a = a s ; P = (3 S ; 7 = 75 • (4) 

There is no s quark in the valence of the proton; hence, we obtain 7y = or 7 = 7 
= 75 (see Table |^). Furthermore, using D^ g = D^ g = D^, leads to the constraint, 
P = 7/4 + a/3. The final simplifying assumption of Du = D-g for all baryons allows 
us to express all antiquark fragmentation functions in terms of 7 alone: 

7 = 7; 

a = 0.757 ; (5) 
P = O.57; 

The model for octet baryons therefore has three unknown functions ay(xE, Q 2 ), 
Pv{ x e,Q 2 ), l(%E,Q 2 ) and an unknown parameter, A, that characterises SU(3) sym- 
metry breaking. We now attempt to evaluate these by suitable comparison with 
data. 

We expect valence fragmentation functions (leading quark fragmentation) to dom- 
inate in the large xe region and sea fragmentation functions to dominate at small 
Xe- We use the large xe data to fix the value of A. Since S~ has two s quarks in its 
valence, its s- valence fragmentation function is suppresesed by a factor of A compared 
to p, A and T, ± '°. Therefore, we expect S~ production cross sections to be smaller 
than for other baryons even in the large Xe range, where only the valence contribution 
survives. Indeed, as can be seen from Fig. [I], data || show that the large xe cross 
sections are similar for all octet baryons, p, A , S 1 * 1 ' , within errorbars, while H~ data 
is smaller in that region. Furthermore, assuming that u quark fragmentation domi- 
nates the large xe proton data, we see from Table || that the ratio of S~ to proton 
production rates is just A (the expressions for the two are otherwise the same). Using 
the data at y/s = 91.2 GeV, we find that A = 0.07. This result can be corrected due 
to the fact that D p d is not small at large xe', however, we shall see that this value for 
A gives good agreement with data. 

Since A turns out to be quite small, the production rate of strange baryons is 
dominated by s-quark fragmentation. This means that the A rate is sensitive to ay 
while the £ baryons give information on (5y (see Table 0). On the other hand, p 
and S depend on the combination {ay + fty); it is therefore possible to separately 
determine ay and (3y from the data, rather accurately, especially in the valence- 
dominated region of large x, x ~ 0.1. The (SU(3) symmetric) sea is the same for all 
members of the octet, upto overall powers of lambda. 

Note that the measured hadron spectrum is an inclusive one. We take this into ac- 
count, especially for the case of p and A, by defining inclusive fragmentation functions 
in terms of the exclusive ones used so far: 

^inclusive = ^exclusive + 0.52S + + 0.64A ; 
Aindusive = ^exclusive + 1-0S° + 1.0^ + 1.0^,° . 

The multiplying fractions indicate the branching fractions into p and A of the various 
baryons. Here S* and A c decays to the various baryons have been ignored since they 
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are very small. A c data is between 15 to 100 times smaller than A and proton data in 
the overlapping xe range of about 0.3-0.8 ||14|| . The S ± and S~ data is considered to 
be purely exclusive for this reason. The energy (xe) of the daughter-baryon is taken 
to be the same as that of the parent because of the small difference in masses of p, 
A, S and S. 

For the case of the 91.2 GeV Z exchange data, ay, (3y and 7 were evaluated at 
different xe values using parametrisations to the p, A and £ data f|, || and the value 
of A as already estimated; these were then used to predict the cross section for E~. 
We find that the data is fitted very well at all xe as is shown in Fig. [I], where the 
resulting fits to ay, (3y and 7 are also shown. Note that a change in A can alter the 
overall normalisation but not the shape of the distribution, which is a prediction of 
this model. 



4.2 Meson Fragmentation 

The pseudo-scalar meson-octet is a self conjugate octet i.e., the same octet contains 
mesons as well as their antiparticles. The mesons_and their antiparticles are related 
by charge conjugation. This means that = D^. As a consequence of this fact we 
immediately see that a, [3 and 7 are not independent of a, (3 and 7 (see Table Q). Of 
these six quantities, only three are independent. We choose these to be a, (3 and 7. 
Due to mixing with the singlet sector, we do not consider the 77 meson here. 

Just as in the case of S, here we predict equal rates for (tt + + 7r~)/2 and ir°, due 
to isospin invariance. This is also borne out by data ||. We therefore consider only 
the combinations = (m + + m~) for m = tt,K, and (K° + K°) = 2K®. 

We make the usual separation into valence and sea fragmentation functions. As 
before, we reduce the number of unknown functions through various symmetry con- 
siderations. We assume that the sea is SU(3) flavour symmetric, so that = 
and so on. Using this, we have (3 = 7/2 and all sea fragmentation functions are equal 
to S = 2 7 . 

Thus, all valence fragmentation functions can be expressed in terms of the func- 
tion, V, where V is given, for example, by the difference (D u — Z\) in tc + , as 

V = a + (3 - | 7 = a - | 7 , 

and all sea fragmentation functions in terms of 7. Now V and 7 can be determined 
by comparison with (7r ± , K ± and K°) data; since the two sets of K data are not very 
different 0, it is not possible to determine a and (3 individually in this case. 

The assumption (made in the baryon case) that the daughter hadron carries away 
the bulk of the energy of the parent, thus contributing to the same xe bin, enabled us 
to simply add the various hadron fragmentation functions to arrive at an "inclusive" 
hadron fragmentation function. Here, since both 7r's and K^s are very much lighter 
than their decay sources (mostly D mesons and baryons), this assumption is no longer 
reasonable. We therefore merely estimate errors arising from the inclusive nature of 
the data by comparing multiplicities rather than xe distributions. 

The bulk of the contamination of the pion sample which is due to K® — > Tin 



decays is estimated to be about 4% from the ratio of the relative multiplicities ||, 15 
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N n /N Ks ~ 10 and a branching fraction ~ 0.35 for the decay for 7r = 7r + , 7r~, 7r°; hence 
we ignore this. K° L does not decay within the detector. However, there is a substantial 
contribution from charm meson feeddown for the K data sample (from the decay of 
all the Ds); this is estimated (from multiplicity data || [L^]) to be about 16% for 
and about 20% for K° mesons. The contamination from baryons is negligible. 

As before, we include SU(3) symmetry breaking effects: the fragmentation proba- 
bility is suppressed by a factor A whenever a strange quark belonging to the valence of 
the meson is produced. Since the fragmenting quark excites a quark pair rather than 
a diquark pair (as in the case of baryons) the value of A here is not related to that for 
the baryonic sector. We can get bounds for A by using the and K ± multiplicities 
= 17.05 ± 0.43 and 2.26 ± 0.18 respectively g [HJ: The total rates for vr ± and K± 
production (their multiplicities) are related to the first moment of V and S, and the 
parameter A; positivity constraints on V and S (since they are probabilities) then 
require that A < 0.14. A tighter constraint on A will be found in the next section, 
when we apply the model to data at different energies. 

We use a typical value of A = 0.08 to fit V and S over the available xe range of 
and K data. These were then used to predict cross sections for tt° and (K° + K°). 
The data is fitted very well by these functions and is shown in Fig. |l|; the fits to V 
and 7 as determined from ir^ and K are also shown here. 

We have been able to explain the production of the entire meson and baryon 
octet by using SU(3) symmetry and the suppression factor A at y/s = 91.2 GeV. 
Encouraged by this success, we investigate whether the model works for other cm. 
energies also. In the next section, we discuss the evolution of these fragmentation 
functions to different energies. 

5 Comparison with data for photon exchange 
5.1 Leading log evolution 

Over the last decade, several experiments, performed over a wide range of cm. en- 
ergies, (10 GeV < y/s < 91.2 GeV), have reported measurements of cross-section of 
hadron production. It is known that the cross-section is a not a constant over the 
range of cm. energies; for instance, this has enabled the extraction of the running 
coupling constant a s . We use leading order (LO) DGLAP evolution equations JTTJ to 
relate the fragmentation functions (and hence the cross-section) at a given energy Q 
to those at a different energy Q . 

Nonsinglet (including the valence functions ay and (3y for baryons and V for 
mesons) and singlet fragmentation functions evolve differently under evolution. The 
singlet (S h (x, t) = u(x, t) + d(x, t) + s(x, t) + u(x, t) + d(x, t) + s(x, £)) mixes with the 
gluon fragmentation function, g(x, t); here we have used x = Xe, t = Q 2 , and = q 
for convenience. Since the gluon is a flavour singlet, we use one gluon fragmentation 
function for all of the mesons and another one for all the baryons. 

The symmetry between the singlet sector of different baryons is broken by A. 
However, all singlet combinations for the other baryons can be expressed in terms of 
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Tj p and strange valence fragmentation functions. We have 



S A = A£ p + (1 - A)4 ; 

S s = A£ p + (1-A)^; (6) 
£ s ~ = A 2 £ p + A(l-A)s _ 



v 



Similarly, we construct only the singlet combination for the pion. All other meson 
singlets can then be recovered from this, and the valence function V . 

The evolution was done in two steps. First of all, for the baryon octet, the gluon 
fragmentation function and 7 (i.e., all the sea quark fragmentation functions) were 
set to zero at the starting scale, and then radiatively generated. This precludes us 
from having to guess a possible starting gluon function, which is very poorly (if at all) 
constrained. Since the sea of all strange baryons is suppressed by a factor of A = 0.07, 
only the proton data (with an unsuppressed sea sector) was poorly fitted. Fits to 
other baryons were still fairly reasonable. This indicates that the fragmentation of, 
say, A or £ is dominated by strange quark fragmentation at almost all xe values. 
Clearly, the function 7 is necessary to fit the proton data. In Ref. |4j], an analysis 
for fragmentation functions for A and p has been done using SU(6) symmetry of the 
baryons. The functions S and T given therein are analogous to the functions a and 
p. However, there is no function analogous to 7, which (being the sea contribution) 
is required to explain the small x data. 

We then used a small, non-zero starting sea as well as (a common) gluon frag- 
mentation function to improve the fits, the results of which are shown in Figs. |2] 
and || However, we emphasise that the evolved fragmentation functions are not very 
sensitive to the choice of the gluon fragmentation function, which is therefore not 
well-determined in our model. 

Each of the fragmentation functions was parametrised at a starting value of Q 2 = 2 
GeV 2 and evolved up to the final Q 2 (= s) value of the data. The input functions 
Fi(x) = ay,/3y, 7 (for the case of baryons) and F^x) = V, 7 (for the case of mesons) 
were parametrised as 

Fi(x) =Oi(l- x) bt (x Cl )(l + diX + etx 2 ) . (7) 

The parameters a, b, c, d, e for different input fragmentation functions are given in 
Table |. 

We tuned the starting parameters to yield a good fit to the 90 GeV hadroproduc- 
tion data which is essentially via Z-exchange @, ||. We then used the same set to 
predict the rates for a lower Q 2 value which is dominated by photon exchange. The 
resulting fits on the Z pole (-y/s = 91.2 GeV) and a fit to the available baryon data 
sample Jl0| at \/s = 34 GeV for A = 0.07 are shown in Figs. 0a and ||]a. In the meson 



sector, we find that A is constrained to lie between 0.04-0.12, with A = 0.08 giving 
the best fit to the data. The overall shape of meson data is very well realised at either 
energy, as can be seen from Figs. 0b and 0b. The model parameters yield a reasonable 
fit to all baryon and meson data at these two different energies. The meson data is 
better fitted than the baryon data at both the energies. The discrepancy in overall 
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normalisation could be due to the inclusive nature of the measurement (especially 
acute in the case of p and A) and possible energy dependence of the suppression fac- 
tor A. However, we emphasise that our model is fairly simple; its biggest advantage 
is that it predicts the production rates of several mesons and baryons with relatively 
few inputs. 

Recently, the total inclusive charged hadron cross section has been measured at 
LEP at 161 GeV [16|. We know from the multiplicity data at the Z° pole that 
81% (91%) of the charged particle inclusive cross section is from pions (pions plus 
kaons). Specifically, the total charged particle multiplicity at y/s = 91.2 GeV is 
21.4 ± 0.02 ± 0.43 |0|, of which 17.05 ± 0.43 are vr ± and 2.26 ± 0.01 ± 0.16 ± 0.09 
are mesons. We therefore compare the charged particle spectrum at 161 GeV 
(with multiplicity 24.46 ±0.45 ±0.44) with our predictions for and (tt^ + K ± ); we 
expect the latter should saturate the data to within 10%. Our model shows excellent 
agreement with data, as can be seen from Fig. [|. 

We remark that the charge factors c q /J2 q c q f° r quarks q — u,d, s, are very different 
for pure Z° and photon exchange. For instance, the charge factor for an s quark is 
1/6 at 34 GeV and 13/36 at 91.2 GeV, more than a factor of two larger. On the 
other hand, that for the u quark is almost a factor of two smaller. This means that 
the photon exchange data is more sensitive to u quark fragmentation than the Z° 
data. That the model predictions for strange hadrons such as A and K (where the u 
contribution is suppressed by a factor of A) are systematically smaller than data may 
therefore mean that D u is actually larger than the model prediction, thus indicating 
that a single strangeness suppression factor A may not suffice. In other words, our 
simple model may not completely account for all SU(3) breaking effects. In this 
context, it would be interesting to obtain data on the £ baryon at a different energy 
and check whether this trend is visible there as well. 

Finally, the data (specially for mesons, p and A) show a decreasing trend at low 
x. The usual DGLAP evolution [|11|] cannot account for such a trend since the pole in 
the splitting function P gg always drives the gluon, and hence the sea, to larger values 
at small x. In 1988, Dokshitzer, Khoze, and Troyan proposed a model wherein 
this dip could be accounted for by including gluon coherence effects. The resulting 
modified leading log approximation (MLLA) then gives a gaussian distribution for 
the singlet fragmentation functions. In the next section, we discuss singlet evolution 
using MLLA and look for improved fits to the low x data. 



5.2 Modified Leading log evolution 



The main result of the MLLA evolution fll2| , [17j is that the low-x singlet fragmentation 
functions have a gaussian form in the variable log(x p ): 



3*p-D(xp) Q) 



N(Q) 



2vra(Q) 



cxp 



-[log(x p )-log(x )] 2 /[2^(Q)] 



(8) 



where N(Q) is the total multiplicity, o is the width of the gaussian and xq is the posi- 
tion of the peak of the gaussian. The Q 2 dependence of N, a and xq are computable 
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for total inclusive hadrons within this approach. They are given as an expansion in 
terms of the scale parameter, Y = log(Q/A), A = 200 MeV: 



N(Q) oc Y- B/2+1/i expJl6N c Y/b 



a' 



,2 _ \r1 



Y'/(3z) ; (9) 



log(l/x ) = Y 



1/2 + Jc/Y-c/Y 



where N c and nj are the number of colours and flavours, which determine the con- 
stants, 

a = UN e /3 + 2n / /(3iV c 2 ) ; b = ( 11N C - 2n/) /3 ; 

B = a/b; z = y/(16N c Y)/b ; 

c = (11/48) [l + (2n / )/(lliV c 3 )] 2 / [1-2^/(11^)]. 

The total multiplicities (at 91.2 GeV, for instance) [§, || [L5| can be used to fix the pro- 
portionality constant for N(Q); the individual particle multiplicities then determine 
N h (q), the multiplicity of the specific hadron, h. The values of a and xq are in good 
agreement with inclusive data fL2] ; however, we are here interested in semi-inclusive 
spectra. In general, the peak shifts to smaller x (here meaning x p ) values for heavier 
hadrons. Also, the semi-inclusive widths are naturally smaller than the total inclusive 
ones. We therefore parametrise the corresponding semi-inclusive parameters as 



x 



o 



ol, 



Cfx 



where Xq is the peak position and the width of the data for hadron h. Here C± 
and C*2 are Q 2 independent constants which we fit to the 91.2 GeV data. They are 
given in Table |]. These are then used to determine the rates at lower energies. The 
resulting fits are again quite good and shown in Figs. [| and || for 90 and 34 GeV 
respectively. Note that the cross section has been plotted as a function of x v and not 
xe here. 

Note that the MLLA is a fit to the singlet fragmentation functions alone; therefore 
comparison should be made with data for x p ~ 0.1, where the valence contribution is 
expected to be small. In the case of Z° exchange, this is also a good fit to the entire 
data. This is because at 91.2 GeV, the cross section is dominated by the singlet term, 
as can be seen from writing Eq. [3] explicitly: 

1 da h ' z _ l2Z h - 1.5£>3 - 0.5L> 8 

otot dx 36 

We see that the singlet contribution is about 10 times larger than either of the octet 
contributions. We therefore expect the MLLA approach to yield sensible fits to the 
data at this energy. In the case of photon exchange data, the singlet contribution is 
still large: 

1 da fe ' 7 _ 2E h + 1.5£>3 + 0-5£>8 

a tot dx 6 
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but the -D3 contribution is not small. Hence MLLA may not be a very good description 
of the data at smaller energies, especially at larger x. However, for the case of A and 
E°, D3 = so that the MLLA singlet term still saturates the event rate to a good 
approximation. 

6 Discussion and Conclusions 

We have proposed a simple model for quark fragmentation into an octet baryon 
or a pseudoscalar meson, using SU(3) symmetry of quarks and octet hadrons. All 
quark fragmentation functions have been described in terms of three SU(3) symmetric 
functions a(x, Q 2 ), fi(x,Q 2 ) and j(x, Q 2 ) and an SU(3) breaking parameter A. The 
antiquark fragmentation functions are correspondingly described by a, (3 and 7. There 
are 3 quark (plus 3 antiquark) fragmentation functions corresponding to a given 
hadron; hence a given hadron octet would involve a total of (24 x 2) fragmentation 
functions. All these are described in our model by just 6 functions, leading to a 
very simple model, but with strong predictive power. Leading log evolution of these 
fragmentation functions has been used to compare the model predictions with data. 
We find that it is possible to fit the model parameters in such a way as to get a good 
agreement with the ^-dependence of all octet baryons and mesons, at two different 
sample energies (corresponding to Z° and photon exchange) over most of the x range 
of available data. These fits were then used to determine the inclusive cross section 
at 161 GeV where both photon and Z° exchange are involved. There was good 
agreement with data here as well. We have used both DGLAP evolution as well as the 
modified leading log approximation (MLLA) to evolve the fragmentation functions; 
the latter has especially been used to explain the decrease in the hadroproduction 
rates at small-x seen in the data for many hadrons. We have used a small non-zero 
input gluon distribution (as shown in Table |J); however, very little sensitivity to the 
gluon fragmentation function is seen and this is therefore not well-determined in our 
analysis. 

The model realises the shape of the ^-distribution of all available data on octet 
mesons and baryons very well. It does not describe the A and K data at 34 GeV 
very accurately; however, it is able to give a good agreement with even this data to 
within 2a. All other baryon and meson data are fitted very well. However, we note 
that it is possible to get good fits at all Q 2 for each hadron individually. The SU(3) 
symmetry constraint relating the different hadron fragmentation functions worsens 
the fit in some cases; this reflects the simplicity of our model, which incorporates 
SU(3) symmetry breaking effects in a very simple way. The goodness-of-fit from 
the model therefore also indicates the extent to which this symmetry breaking is a 
universal phenomenon, independent of the type of quark or diquark that is produced 

B0.0I- 

The parameter A takes into account the difference in masses of the strange and non- 
strange quarks and suppresses non-strange quark fragmentation into strange hadrons. 
This parameter is similar to the suppression factor of the Lund Monte Carlo [|J, 
however, it is determined by means of a simple comparison of data of stange and 
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non-strange hadrons. The Lund model uses string fragmentation and has a much 
larger suppression for the case of baryons (suppression factor = 0.06) as compared to 
the suppression factor of 0.2-0.3 for mesons. Our model has very similar values for 
the suppression factors for the two cases (A = 0.07 for baryons, 0.08 for mesons). 

Another approach Q uses an SU(6) analysis of fragmentation functions using a 
quark and diquark model. Our SU(3) symmetric functions a and (3 are analogous to 
the SU(6) symmetric functions S(z) and T(z) defined therein. The function 7 (not 
included in their model) describes sea fragmentation, for instance, s fragmenting to 
a proton. We find that 7 is large in the small x region, so that its contribution is 
significant and cannot be ignored. 

We find that the strange quark fragmentation dominates strange hadroproduction 
over almost the entire x range. This is especially true for A, which has recently been 
of much interest Jl|, 18]. It is possible to extend the model to include spin-dependent 
fragmentation functions; the unpolarised result then indicates that polarised A frag- 
mentation will be dominated by its strange fragmentation function, which can then 
be readily parametrised and studied. 

Finally, our results suggest that there is indeed an underlying symmetry among 
the baryons and mesons in an octet, which can be tested further by extending the 
model to decuplet baryons and other hadrons. 
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Appendix A 

We briefly detail the calculations leading to the results in Table ^| for the quark 
fragmentation functions in terms of a, (3 and 7. 

Let qi be a quark triplet, B\ the hadron octet, k = 1, 2, 3. 

Case 1 A is a triplet, Xf. Then the invariant amplitude for the process q — > H + 
X is qiHjXi , where Aj and are normalised. Here H\ are the elements of the 
meson/baryon matrix (see Table [l]). Thus, the rate for u — > p + A is a\uBlX 3 \ 2 
which is equal to a. Similarly, the rate for u — > A + A is a\uB\X l \ 2 which is equal 
to a/6 and so on. 

Case 2 A is a sixplet, A^-: Now Ay is symmetric in % and j and is expressed 
in terms of triplets as (q%qj + qjqi)/V2, where each q is normalised. The invariant 
amplitude is e imj qi H*X km . Thus, d -> p + X as (3\y/2\ 2 and so on. 
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Case 3 X is a fifteenplet, Xj k : Then X( k is symmetric in j, k and is antisymmetric 
in In terms of triplets, the normalised X can be re-expressed as 

q j q k qi + q k q j qi ~ ~^l{q 1 q k qi + q k q l qi) - ^S k {q j q l qi + q l q j qi ) , 

where each of the q^s is normalised. 

The invariant amplitude is X^Hjqk. In this case, we shall have to take the 
interference terms for the diagonal elements of the meson/baryon matrix also into 
account. Note that X\ 3 = (sum over i is implied). Thus, the rate for u —> A + X is 
^(l/V^Xl 1 + (l/VQ)X? - (2/v / 6)X| 1 | 2 which is equal to 7 |(3/ V / 6)(X 1 11 + Xf )\ 2 . 
On evaluating this expression, we find that this is equal to 9/87. The other rates can 
also be found in a similar manner. The final results are given in Table |2|. 
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Table 1: (a) Members of the meson octet and (b) Members of the baryon octet. 
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Table 2: Quark fragmentation functions into members of the baryon and meson octet 
in terms of the SU(3) functions, a, /3 and 7. 
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Pv 


7 
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a 
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Table 3: Input values at Q 2 = 2 GeV 2 for the valence and singlet fragmentation 
functions for the (a) meson and (b) baryon octet. 
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0.92 
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0.75 


c 2 


0.4 


0.4 


0.4 


0.4 


c 2 


0.59 


0.59 


0.425 


0.425 



Table 4: Multiplying constant factors for description of singlet fragmentation func- 
tions for mesons and baryons in the MLLA approach. 
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Figure 1: The figure on the left shows the baryon fragmentation probabilities acy, ftv 
and 7, fitted using data on p, A and E at 90 GeV || |9| as a function of x, and the 
prediction for the S baryon using these, along with the data used. The figure on the 
right shows the mesonic probabilities, V and 7, and the predictions for mesons; the 
fits shown are to data corresponding to 7r ± , 7r°, K ± and K° in decreasing order of 
magnitude. 
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Figure 2: (a) The figure shows the model fits to the 90 GeV data for baryons using 
DGLAP evolution. 
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Figure 2: (b) The figure shows the model fits to the 90 GeV data for mesons using 
DGLAP evolution. 



19 




Figure 3: (a) The figure shows the model fits to the 34 GeV data for baryons using 
DGLAP evolution. 




Figure 3: (b) The figure shows the model fits to the 34 GeV data for mesons using 
DGLAP evolution. 
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Figure 5: (a) The figure shows the model fits to the 90 GeV data for baryons using 
MLLA evolution. Note that we have used x = x p in order to clearly exhibit the 
small-x data which is of interest here. 
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Figure 5: (b) The figure shows the model fits to the 90 GeV data for mesons using 
MLLA evolution. Note that we have used x = x p in order to clearly exhibit the 
small-x data which is of interest here. 
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Figure 6: (a) The figure shows the model fits to the 34 GeV data for baryons using 
MLLA evolution. Note that we have used x = x p in order to clearly exhibit the 
small-x data which is of interest here. 




Figure 6: (b) The figure shows the model fits to the 34 GeV data for mesons using 
MLLA evolution. Note that we have used x = x p in order to clearly exhibit the 
small-rc data which is of interest here. 
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